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ABSTRACT 

We present radial velocities for 38 low-mass candidate members of the cr Orionis 
young group. We have measured their radial velocities by cross-correlation of high 
resolution (R« 6000) AF2/WYFF0S spectra of the gravity sensitive Na i doublet at 
8183, 8195 A. The total sample contained 117 objects of which 54 have sufficient signal- 
to-noise to detect Na i at an equivalent width of 3 A, however we only detect Na i in 38 
of these. This implies that very low-mass members of this young group display weaker 
Nal absorption than similarly aged objects in the Upper Scorpius OB association. 
We develop a technique to assess membership using radial velocities with a range of 
uncertainties that does not bias the selection when large uncertainties are present. 
The resulting membership probabilities are used to assess the issue of exclusion in 
photometric selections, and we find that very few members are likely to be excluded 
by such techniques. We also assess the level of contamination in the expected pre- 
main sequence region of colour-magnitude space brighter than I = 17. Wc find that 
contamination by non-members in the expected PMS region of the colour-magnitude 
diagram is small. We conclude that although radial velocity alone is insufficient to 
confirm membership, high signal-to-noise observations of the Nai doublet provide the 
opportunity to use the strength of Nai absorption in concert with radial velocities 
to asses membership down to the lowest masses, where Lithium absorption no longer 
distinguishes youth. 

Key words: techniques: radial velocities - stars: pre-main-sequence - stars: low- 
mass, brown dwarfs - techniques: spectroscopic - open clusters and associations: a 
Orionis - 



Studies of the low-mass and sub-stellar initial mass func- 
tion (IMF) provide the opportunity to constrain theories of 
star formation and investigate environmental influences on 
the final results of the formation process. Whilst the high 
and intermediate mass regions of the IMF have been well 
studied and found to be essentially universal in the field 
and in clusters (e.g. lHillenbrand 2003t ICarmanv et al. 1982L 
and references therein), this is not the case for low and sub- 
stellar masses. Although there have been a number of stud- 
ies of the low-mass IMF in the field (e.g. iReid et al. 19991 : 



1 INTRODUCTION 



iKroupa 200lh . due to their inherent faintness only the most 
nearby of field T dwarfs are detected, leading to poor statis- 
tics at these masses. Searches for brown dwarfs in interme- 
diate age (10* yrs) clusters such as the Pleiades have proved 
fruitful (e.g. ^oraux ot al. 2003; Basri 2000, and references 
therein). Unfortunately, to derive an IMF from the mass 
function of an intermediate age cluster requires corrections 
for the effects of stellar and dynamical evolution over the 
age of the cluster. 



Young clusters and associations are the ideal places 
to study this region of the IMF. This is because many of 
the problems associated with measuring the IMF at the 
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lowest masses disappear in such young regions. Firstly no 
correction need be made for the effects of stellar or dy- 
namical evolution, the regions are sufficiently young that 
the measured mass function is the IMF. Also, the ob- 
jects at substellar masses are much easier to observe at 
ages of less than 10 Myrs since they still shine brightly 
as gravitational potential is released. In addition, a sup- 
posed advantage of clusters and associations with ages of 
less than 2 Myrs is the lack of contaminating sources due 
to the fact that such young clusters are often still associ- 
ated with their natal molecular cloud, which will tend to 
obscure the background to high extinctions. The Orion Neb- 
ula Cluster (ON C) has been extensively studied for just 
these reasons Ce.g.lMuench et al. 20o1lLuhman et al. 2000l: 
iLucas fc Roche 200ffl . However, inspection of the methods 
which must be employed for such studies reveals that asso- 
ciation with a molecular cloud brings with it its own, un- 
avoidable, problems. 

Since there is often a spread of ages in young star form- 
ing regions, it is preferable to use photometric magnitudes 
and spectroscopically derived temperatures in conjunction 
with pre-main sequence (PMS) models to determine the age 
and mass for each object in the sample. Unfortunately due to 
the high extinctions present along sightlines toward young, 
partially embedded, clusters it is generally not possible to 
obtain spectroscopy for the lowest mass objects. Since there 
is also normally a large range of reddening toward objects 
within such clusters it is not possible to determine the age 
and mass for each object, and a distribution of ages and 
reddening must be assu med to obtain the IM F from the lu- 
minosity function (LF) (iMuench et al. 20021 ). Another effect 
of the high extinction is that LFs must be constructed in the 
K- and L-band to avoid biasing the sample to the least ex- 
tincted objects and to ensure a good level of completeness. 
However, it has been observed in the young clusters such 
as the ONC that up to 85% of the objects have circumstel- 
lar discs, which are responsible for infra-red excesses (e.g. 
iLada et al. 2000l) . making mass functions derived from K- 
and L-band LFs unreliable. These problems are avoided if 
star forming regions between 2 and 10 Myrs old are used for 
IMF studies. At these ages the clusters have often emerged 
fully from their natal cloud, presumably due to the UV flux 
from the massive stars evaporating the molecular gas. The 
lower extinction allows the luminosity function to be con- 
structed in the I-band, which is less affected by circumstellar 
discs than the K-band, in addition to making spectroscopy 
possible down to very low masses. Since we need such re- 
gions to be as nearby as possible, to allow studies of the 
lowest mass objects, the number of prospective laboratories 
is very small. The most promising candidate is the o Orio- 
nis (cr Ori) young group, a constituent of the Orion OBlb 
associ ation, w hich has age estimate s in the ran ge 1.7-5 
Myrs JWarren fc Hcsscr 1978; Bro wn et al. 1994) . The low 
extinction, E{B — V) = 0.05 llLee 196'8[) . and small distance, 
D « 350pc jPerrvman et al. 1997^ . make it ideal for studies 
of the low mass IMF. 

iBeiar et al. f200ll) have carried out an investigation of 
the substellar IMF in the a Ori region. They used a sam- 
ple of 9 spectroscopically confirmed members to define the 
locus of the cluster sequence between M6 and L4 spec- 
tral types. They then used the location of this sequence 
in colour-magnitude (C-M) space to identify cluster mem- 



bers, and construct an l-band L F. Using several d ifferent 
models to construct their IMF, iBeiar et al. f200J) found 
a mass spectrum (dN/dm) oc m~°' with an exponent of 
a =0.8±0.4 between 0.2 and O.OIM©, i.e. the number of 
objects per mass bin continues to increase with decreas- 
ing mass down to below the deuterium burning limit and 
into the realm of planetary mass objects. This corresponds 
to an IMF ^(m) oc M"' with an exponent 7 = -0.2. The 
low extinction, however, means that background contami- 
nation ca nnot be ruled out, even in the locu s of the cluster 
sequence. iBarrado v Navascues et al. (200 J) obtained spec- 
tral types for roughly one quarter of the photometrically 
selected candidates, confirming a high proportion of them 
as cool objects. Whilst such a method is effective at remov- 
ing reddened background objects from a sample, it is not 
robust against interloping field dwarfs. 

When the number of objects in each mass bin is so 
small, however, a small number of contaminants can have 
a significant effec t on the derived mass function. Whilst 
iBeiar et al. f200J) are careful with their photometric selec- 
tion to avoid contamination, this could, in itself, lead to 
poorly understood completeness. Without empirically deter- 
mining how far into the expected background region of the 
colour-magnitude diagram (CMD) bona fide members are 
still found, any photometric selection which minimises con- 
tamination may well exclude a significant fraction of member 
objects. 

The likelihood of age spread in such a young cl uster 
makes this consideration crucial. iKenvon et al. f2004) ob- 
tained fibre spectroscopy for a sample of over 70 candidate 
members, drawn from a RI catalogue, with masses between 
0.055 and O.SMq (they assume a cluster age of 5 Myrs). 
They found that contamination from background objects 
was limited for the reddest region of each magnitude bin 
in the I/R-I CMD, but the PMS objects were poorly cor- 
related with isochrones in this colour. They found bette r 
correlation in I/I-J, the colours used bv lBeiar et al. f200ll) . 
but the contamination was found to be considerably worse, 
becoming more prevalent at the faintest magnitudes. Addi- 
tionally, the CMDs of Konvon ct al. (2004) show that any 
strict cut along isochrones in either set of colours would ex- 
clude a significant number of members. 

It is clear that the only way to obtain a reliable IMF for 
this region is to use a spectroscopically confirmed sample of 
memb ers, drawn from a br oad region of colour magnitude 
space. iKenvon et al. (2004^ used the presence of Li I ab- 
sorption at 6708 A to confirm youth and thus membership 
for objects in their sample. This technique, however, is lim- 
ited to objects with M > 0.065Mq . Since objects below this 
mass never attain high enough core temperatures to burn 
Lithium, old foreground T-dwarfs cannot be distinguished 
from bona fide members based on Li absorption alone. 

If a reliable IMF is to be derived for the lowest masses, 
another method of distinguishing members must be found. It 
is with this in mind that we have obtained fibre spectroscopy 
of the Nai doublet at 8183, 8195 A for a photometrically se- 
lected sample of objects in the direction of the a Ori cluster. 
Using this doublet has several benefits. Cool objects such as 
low-mass stars and brown dwarfs are brightest in the NIR 
and IR regions, where night sky emission can severely affect 
observations. The Nai doublet lies at a wavelength which is 
relatively unaffected by bright sky emission lines. Also the 
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Nal doublet is a relatively strong spectral feature in cool 
stars, and so should be easy to observe a moderate signal- 
to-noise. 

In this paper we have two aims. Firstly, to investigate 
the use of radial velocity, as measured from the Nal doublet, 
as a membership diagnostic. Whilst Li absorption would be 
an effective diagnostic for much of our sample, this work is 
concerned with developing a spectroscopic technique for use 
at magnitudes and colours where foreground contaminants 
may lie below the Li burning limit. Our second aim is to 
measure the level of contamination from non-members in the 
expected PMS region of the CMD and to determine if many 
bona fide members are excluded by photometric selection 
techniques. 

The paper will be laid out as follows. In Section |21 we 
will discuss the rationale behind our target selection and 
our observations. In Section |3 we will describe our data re- 
duction technique and our method for sky subtraction. In 
Section 0] we will describe our sample selection, the cross 
correlation of spectra to obtain radial velocities and the cal- 
culation of membership probabilities. Section|5]will contain 
a discussion of the results. A summary of our conclusions 
will be given in Section |S| 



2 OPTICAL SPECTROSCOPY USING 
AF2/WYFFOS 

2.1 Target Selection 

We obtained spectra for 117 objects using the 
AF2/WYFFOS fibre spectrograph on the William Herschel 
Telescope at the Observatorio del Roque de Los Muchachos, 
La Palma. The targets were select ed from the photo - 
metric RI catalogue described by iKenvon et al. (200411 
and were chosen to compliment the sample for which 
spectroscopic obser vations have already been obtained 
jKenvon et al. 20041) . Since many of the reddest objects 
have already been observed, particularly at the brighter 
magnitudes, we were able to cut deep into the region which 
is thought to be background contamination, see Figure 
This method of target selection has served a dual purpose: 
1) it allows us to investigate the possibility of contamination 
in the PMS region of colour magnitude space; 2) it allows 
us to investigate the number of member objects which 
might be lost in more conservatively drawn samples. The 
distribution of the targets, in relation to a Ori, and the 
WYFFOS/WHT FoV is shown in Figure H 

2.2 Observations 

We observed our targets using the small fibre (1.6") setup on 
the nights of the 3'^'' and 4"* January 2003. Conditions were 
dark, and cloud was only present during the first half of the 
second night. We observed in 3rd order echelle mode, with a 
central wavelength of 8300 A. This gives a dispersion of 0.57 
A/pixel and we achieved a resolution of approximately 1.4 A. 
We elected to observe the Nai doublet at high resolution so 
as to avoid broadening the bright sky emission lines and 
obscuring the doublet. 

We achieved a total exposure time of 9.5 hours for the 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 
R-I 



Figure 1. The CMD for the optical RI catalogue from which 
targets were selected iKenvon et al. 2004'). Targets for this sur- 
vey are sh own as circles. The dotted line follow s a NextG en 5 Myr 
isochrone iChabrier fc Baraffe 19971: ISaraffe et al. 20021) . We de- 
fine the expected background region to be blueward of the solid 
line, and the PMS region to be redward of it (see section 5). 
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Figure 2. The locations of our targets and WHT FoVs. cr Ori is 
shown as star, whilst our targets are shown as circles. The WHT 
FoVs are shown as large circles. 
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faintest targets. This time was split between 4 configura- 
tions to maximise the number of objects we could observe 
(see Table 0. The brightest targets were only observed 
in the shortest set-ups, with fainter ones observed using a 
combination of set-ups to achieve appropriate integration 
times. Each set of exposures was subdivided into Sections 
of 900s or 1800s, depending on the total exposure required 
for the set-up. Between exposures of the target objects, 
the telescope was moved such that o Ori lay over a fibre, 
and a spectrum of a Ori was obtained for the purposes 
of telluric correction. Copper-Neon arc frames were also 
obtained for both positions of the telescope at these times. 
To maximise our exposure ti mes we did not observed offset 
sky fields, as is advised by IWvse &: Gilmore ('199 2'). for 
estimating the relative throughput of the fibres prior to 
sky subtraction. However, the relative throughput of the 
fibres in the small fibre bundle is considerably more homo- 
geneous than it had been for the older large fibre bundle (see 

:/ /www. ing.iac.es/Astronomy/instruments/af2/index. html I, 
and we have an alternative method for achieving sky sub- 
traction, which will be described later. We also observed 
several stars for use as radial velocity standards, these ob- 
servations are detailed in Table |5| along with the identifier 
of the telluric reference observed. 

Unfortunately, many of the spectra we obtained sufi'ered 
from poor signal-to-noise, and others had no detectable sig- 
nal at all. We have correlated the RI photometric catalogue 
used to select our targets with a 2MASS catalogue for the 
same region and find the RMS in the residuals of the coor- 
dinates to be approximately 0.2" in both axes. As such, we 
rule out errors in the astrometry of the RI catalogue as an 
explanation for our missed targets. We also rule out proper 
motion as a factor. At the target selection stage the 2MASS 
catalogue for our potential targets was correlated against 
the SuperCosmos sky survey catalogue for the same region 
and objects with large coordinate residuals were excluded. 
Whilst source variability, well known in T Tauri stars, may 
have had an infiuence on the signal detected from some of the 
targets, we believe the failure to achieve the desired signal- 
to- noise in so many spectra is the result of two problems with 
the AF2/WYFFOS field of view (FoV). Firstly, the point 
spread function (PSF) is highly variable across the FoV, be- 
coming very smeared out beyond a radius of 25 arcminutes. 
As such, when using the small fibre bundle a large propor- 
tion of the signal from a target star in the outer area of the 
FoV will be lost. This effect explains the poor signal-to-noise 
found in the spectra of some objects. The second problem is 
that the astrometric distortion in the AF2/WYFFOS FoV 
is poorly constrained, meaning that many objects were sim- 
ply missed by their assigned fibre. The method employed for 
selecting our sample from the available spectra is discussed 
in Section f4.1l 



3 DATA REDUCTION 
3.1 Extraction of Spectra 

The spectra were extracted and dispersion corrected using 
the WYFFOS data reduction routi nes within th e RGO pack- 
age within the IRAF environment ijLewis 199(t) . The spectra 
were bias subtracted using the mean bias level in the over- 



Table 1. Summary of the exposure times and sky fibre allocations 
for the 4 configurations used. 



Configuration 


Integration Time (s) 


No. of Sky Fibres 


Longl 


12600 


58 


Long2 


12600 


49 


Medium 


5400 


35 


Short 


3600 


29 



scan region of the Tek6 detector. Flatfielding was achieved 
using Tungsten fiatfields, which were also used for aperture 
detection. The apertures were traced using steps of 5 pixels 
and fitted using a 2'"'' order spline function. Spectra were 
then optimally extracted from target, g Ori and arc frames. 
Dispersion correction for each spectrum was achieved by first 
fitting a 4"* order Chebyshev function to arc lines that were 
manually identified in the spectrum from the middle aper- 
ture of the appropriate arc frame. An automated procedure 
then obtained solutions using the same lines in the remaining 
apertures. The RMS deviation from the dispersion solution 
across the apertures in the arc frames was typically 0.01- 
0.02 A. These solutions were then applied to spectra from 
the same apertures in the target frame. 



3.2 Telluric Correction and Sky Subtraction 

Absorption features due the Earth's atmosphere can influ- 
ence the results of cross-correlation, by biasing the result to- 
wards zero velocity. As such the effects of the atmospheric, 
or telluric, absorption must be corrected for. It is normal 
for an 0-star, whose spectrum should be nearly featureless, 
to be used to as a reference spectrum for telluric correction. 
This can lead to difficulty when the reference star is well sep- 
arated from the targets, as one is required to assume that 
the telluric absorption is uniform in both time and position. 
Fortunately, a Ori itself provides us with an appropriate ref- 
erence spectrum, and is located within our FoVs. As such, 
we need only assume that the telluric absorption is uniform 
over short periods of time. 

Telluric correction was carried out on spectra from each 
sub-exposure using the temporally closest o Ori spectrum. 
Radial velocity standard spectra were corrected using a tel- 
luric reference observed immediately before or after each 
exposure, and chosen to be close-by on the sky, and these 
are listed in Table |5| In each case the object spectrum was 
divided by a rectified version of the telluric reference spec- 
trum, thus removing absorption features originating in the 
Earth's atmosphere. 

The object spectra from the sub-exposures of each con- 
figuration were combined by median stacking prior to sky 
subtraction. We calculated the uncertainties in flux at this 
point as the standard error about the median. Two master 
sky spectra were constructed for each configuration by calcu- 
lating (i) the weighted mean and (ii) the median of median 
stacked spectra from each sub-exposure. The sky spectra 
used for constructing the median sky for each sub-exposure 
were selected on the basis of being free of cosmic rays in the 
region of interest (8100 A < A < 8300 A). Since our target 
region is well seperated from the nebular emission near the 
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Horsehead, the background is flat. As such, use of a master 
sky spectrum for each set-up is appropriate (see below also) . 

Since no offset sky frames were obtained, it was not 
possible to estimate the relative throughput of each fibre 
and use this to scale the master sky spectrum to each fi- 
bre for subtraction. Instead we have iteratively scaled the 
master sky frame until the subtraction is optimised for the 
removal of sk y lines. The meth od we use is identical to that 
used by Beekman et al. f200Cf) for removing the, relatively 
sharp-spiked, secondary star spectrum from that of a cat- 
aclysmic variable (CV). In this case the master sky frame 
is treated as a template secondary star spectrum, and the 
target+sky as the CV spectrum. First we subtract the mas- 
ter sky frame with no scaling and smooth the residual with 
a smoothing length of 50 pixels. The smoothed residual is 
then subtracted from the unsmoothed residual and the 
per pixel calculated. This procedure is then repeated with a 
series of different scale factors. The scale factor that results 
in the lowest value of is the one which is then applied to 
master sky spectrum prior to subtraction from the data. An 
example of the success of this method of sky subtraction is 
shown in Figure |H] This method was carried out using both 
versions of the master sky frames. As can be seen in Figure|21 
sky emission lines lie at similar wavelengths to the Nal dou- 
blet of interest. As such it is crucial that the sky lines are 
not over-subtracted, giving rise to spurious absorption fea- 
tures. To assess the effectiveness of each sky subtraction we 
used the two emission lines just blueward of 8150 A as diag- 
nostic lines. Since these lines were found to be more intense 
than those situated over Nal doublet in the spectra from all 
sky fibres, they are good indicators of over-subtraction. A 
problem would arise if the relative strengths of these lines 
varied significantly across the field of view, or over the pe- 
riod of time that a set-up was observed for. To assess this 
potential problem we have tested the sky subtraction on a 
number of sky fibres from across the field of view and from 
different exposures in a set-up. No cases were found where 
over-subtraction had taken place but was not evidenced by 
absorption in the diagnostic lines. Figure |l]shows a selection 
of these test subtractions to illustrate this point. 

The visual inspection of the residual spectra resulted 
in a preferred version of the sky frames to be used for each 
object in a configuration, and this was forced to be the same 
across all configurations for a given object to ensure a con- 
sistent method. In the vast majority of cases the median 
stacked sky frames performed better. It should be noted that 
since the spectra have not been corrected for the relative 
throughput of the fibres, it is not possible to fiux calibrate 
spectra that have been sky subtracted in this manner. At 
this stage spectra for objects obtained in different configu- 
rations were co-added to arrive at the final spectra. These 
spectra are available on-line via the CDS service or from the 
Cluster collaboration homepage^. 
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Figure 3. An example of the effectiveness of our method for sky 
subtraction. In this case, the median stacked exposures from 3.5 
hours of the 7 hour total exposure time for our target 1.03-342 
are sky subtracted using the method described above. 
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Figure 4. Examples of sky subtracted sky spectra used to check 
validity of using the diagnostic lines near 8150 A 



4 SELECTING MEMBERS BY RADIAL 
VELOCITY 



4.1 Sample Selection 



As discussed in Section 12.11 we did not obtain good data 
for every object in the sample. To avoid biasing our sam- 
ple at this stage we opted to apply a signal-to-noise cut 
to select our sample. We measured the equivalent width of 
the Nal doublet, hereafter EW(Nal), in all the spectra ob- 
tained. We used two continuum bands: 1) 8151-8173 A; 2) 
8233-8255 A. The wavelength range for continuum band 2 
was selected to avoid the TiO band which lies redward of 
8200 A. We integrated EW(Nai) over the range 8180-8202 A. 
Those objects which had a sufficiently small uncertainty to 
allow a 2a detection of the doublet at an equivalent width 
http://www. astro. ex.ac.uk/people/timn/Catalogues/description.html(EW) of 3 A, i.e. a > 1.5 A, were included in the sam- 
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pie. This EW was c hosen as most brown dwarfs detected by 
iMartfn et al. in the Upper Scorpius OB association 

displayed an EW(Nal) above this value. The selected objects 
and their values for EW(Nal) are listed in Tabled Our ob- 
jective measure for the presence of Nal absorption was an 
EW(Nal) measured at a significance of 2cr or greater. Our 
objective measure failed to detect Nai in 11 of the selected 
sample objects. We rule out the presence of Nal in a fur- 
ther 5 objects, despite the objective measure suggesting its 
presence. These false detections were caused in two cases 
(1.03-460, 1.03-612) by the small sky line residuals affecting 
the continuum estimate. One false detection (8.03-396) was 
caused by over-subtracted sky lines within the EW integra- 
tion band. The other two false detections (8.04-50, 8.04-77) 
were caused by unidentified absorption features encroaching 
on the EW integration band. Since no radial velocity can be 
obtained when Nal is not detected, we do not include such 
objects in our cross correlation nor do we calculate member- 
ship probabilities. 

All of the objects with a failed Nal detection are fainter 
than 1=17. The strength of the Nal doublet is gravity de- 
pendant, getting stronger with increasing surface gravity. As 
such it has been used as a discriminato r between field dwarfs 
and giants (e.g. ISchiavon et al. 1997ft an d field dwarfs and 
young brown dwarfs l|Ma^^^^I^Oo3) . Failure to detect 
Nal in the spectrum of an object could indicate that it is 
a giant with a sufficiently low surface gravity to weaken 
Nal beyond our detection threshold. However, since very 
young brown dwarfs also have low surface gravity the same 
might be expected from members of the a Ori group. Sur- 
face gravities are typically higher in young (age < 5 Myrs) 
brown dwarfs, Ioq{q ) ~ 3.5, than in M-giants, log(g) ~ 
iGorlova et al. 2003ft . so we might expect to be able to dis- 
tinguish between extremely low surface gravity giants and 
young brown dwarfs. This will be discussed further in Sec- 
tion |K| 



4.2 Cross Correlation of Spectra to Obtain Radial 
Velocities 

We obtained radial velocities for objects in our sample which 
displayed Nal absorption by cross correlating their spectra 
against those of stars with known velocities and similar spec- 
tral types to our targets. The stars observed as radial veloc- 
ity standards are detailed in table |5| The wavelength range 
used for the cross correlation was restricted to 8175 - 8205 A 
to minimise the effects of noise in the spectra. We restricted 
the range of acceptable correlations to be zblOOkm s~^, and 
as a result velocities were not found for some objects on some 
cross correlations. Uncertainties for each velocity measure- 
ment were estimated by perturbing the value of the flux at 
each data point in the target spectrum by some value drawn 
from a Gaussian based on the uncertainty in the flux. This 
was repeated 100 times and the resulting spectra were cross 
correlated against the standards in the same manner as the 
original data. The uncertainty in velocity was taken to be 
the standard deviation in the distribution of velocities from 
the perturbed spectra. The quality of all cross correlation 
functions were assessed by visual inspection and none were 
found to have any evidence of the multiple peaks that might 



be expected if the sky subtraction introduced spurious fea- 
tures. 

Although the 4 different velocity standards gave similar 
results, the use of GJ3517 gave the most successful corre- 
lations. As such, this is the one we chosen to use for fur- 
ther analysis. The velocities and uncertainties resulting from 
cross correlations against GJ 3517 are given in Table |H] 



4.3 Membership Probabilities 

A histogram of the velocities obtained in the previous Sec- 
tion (Figure 01 shows that there is no simple velocity cut 
that will separate members of the a Ori group from non- 
members. Also the velocities derived in the previous Section 
have a range of uncertainties, due the range of signal-to- 
noise in the different spectra, so a simple na-cut from a 
mean cluster velocity will tend to preferentially select poor 
signal-to-noise objects as members. Instead we can assign a 
probability of membership, Pmem, to each object for which 
we have a velocity. This probability is the product of the 
probability that the object is at an appropriate velocity to 
be considered a member, P„ei, and the probability that an 
object at that velocity is not a radial velocity contaminant, 
(1 - Pcont)- The first question that must be answered, then, is 
what constitutes an appropriate velocity for membership of 
the cluster? We have plotted the cumulative probability dis- 
tribution of the sample by summing Gaussians constructed 
from the velocities and uncertainties for all objects with a 
radial velocity. As can be seen in Figure |S] there is a strong 
peak in this distribution centered at 29.5 km s^^. Since we 
do not resolve the velocity dispersion of the cluster, we can 
take the full width at the level of the background as the 
range of velocities occupied by cluster members. Pvei is then 
simply found as the fraction of the Gaussian derived from an 
object's velocity and error that lies within the cluster range. 
We take the cluster range to be 24-37 km s~^, based on the 
peak's full width at the level of the background. 

The value we find for the group velocity, 29.5 km s~^ 
(and the ra nge of 24-37 km s~^ ), is consistent with that 
measured bv lKenvon et al.' (2004ft of 31.2 km s ^, and with 
the mean of values for members of this group measured by 
iMuzerolle et al. f2003ft of 30.9 km s ^. It is also consistent 
with the mean radial v elocity of a Orionis itself of 29.2 ±2.0 
km s" -1 JWilson 1953ft. This is at odds with the value mea- 
sured bv IZeroater^^sorio et al" r2002ft of 37.3 km s~\ but 
it is not clear why they have measured such a high velocity 
for this young group. 

The probability that an object at the velocity of the 
cluster is a radial velocity contaminant can be taken as the 
ratio of the probability in Figure |S| integrated over a ve- 
locity interval equal in size to the cluster range (but lying 
in the non-cluster region) to the probability integrated over 
the cluster peak in Figure |S| We find Pcont = 0.2 for this 
sample. We emphasise that this probability applies to the 
entire sample, i.e. over the entire C-M space examined. The 
value of Pcont will change if, for example, one limits ones 
attention to the "PMS" region of the CMD (see below). As 
such we only list values of Pvei in table |2l so as to avoid 
misinterpretation of the results. 
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Table 2. Details of stars observed as radial velocity standards and details of their observations. Radial velocities listed are barycentric. 
Equivalent widths for Nal are measured with continuum bands at (8151-8173 A) and (8233-8255 A), and integrated over the range 
(8180-8202 A). 



Star 


a (J2000) 


(5(J2000) 


Sp. Type 


I 


Integration time (s) 


Vrad / km s 1 


EW(Nai) / A 


O'EW 


Telluric reference 


WX UMa 


11 05 30.31 


4-43 31 16.6 


M6 




200 


68^ 


8.28 


0.10 


SAO 43460 


GJ412A 


11 05 28.58 


-1-43 31 36.4 


M0.5 




10 


68i 


2.18 


0.09 


SAO 43460 


GJ3517 


08 53 36.11 


-03 29 32.4 


M9 


14.44 


400 


92 


5.04 


0.38 


SAO 136434 


LP213-67 


10 47 12.65 


-1-40 26 43.7 


M6.5 


13.17 


600 


53 


6.55 


0.21 


SAO 62257 



References for radial velocities: lI lDelfosse et al. (199^) : 2') lTinnev fc Reid (1993) : 3 lReid et al. (20021) . 
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Figure 5. A histogram of the velocities measured from our sam- 
ple. 



5 DISCUSSION 

In the following discussion we refer to the region redward of 
the solid line in Figure |7| as the expected PMS region. The 
region that lies blueward of this line is referred to as the 
background region. This locus was selected as the dividing 
line as it corresponds to a liberal photometric selection, and 
so provides an interesting case for investigating exclusion of 
members and contamination by non-members. 



5.1 Are there members outside the "PMS" 
region? 

The simplest way to answer this question is to examine the 
cumulative probability distribution of the velocities for the 
objects in the "background" region of the CMD. We have 
plotted this distribution with a dotted line in Figure |S| It 
is clear from this plot that there is no over-density in the 
probability distribution at the cluster velocity. This implies 
that there is not a significant number of members in the 
background region of the CMD, and hence photometric se- 
lection techniques are not excluding members of the cr Ori 
young group. This is the principal conclusion of this work. 



Figure 6. A plot of the cumulative probability distribution of 
our sample. The strongest peak is centered at 29.5 km s~^. The 
dotted line traces the cumulative probability distribution of the 
objects in the "background" region of the CMD (see Figure ItI. 



5.2 Contamination within the "PMS" region 

As can be seen in Figure |7| at magnitudes brighter than 
about I = 17 we find that all of the objects in the ex- 
pected PMS region have P„ei > 80% except one, which 
has Pyei = 77.3%. This is consistent with the result of 
iKenvon et al. (2004^ that essentially all objects in the PMS 
region of the CMD are members of the a Ori young group. 

Fainter than I = 17 the situation is somewhat differ- 
ent. We only detect Nal in 13 out of 29 objects in this part 
of our sample. The failure to detect Nal in so many ob- 
jects in this part of the CMD makes a discussion of the 
contamination problems in this region impossible. However, 
this failure is itself worthy of some discussion. As already 
mentioned, the signal-to-noise thresh old used to select ou r 
sample was chosen by comparison to iMartm et al. (20041) ■ 
who find EW(Nal) > 3 A for all but 2 of their members 
of the 5 Myr old Upper Scorpius association. As such we 
expected that the majority of members of the similar aged 
(J Ori group might share this property. By comparing the 
EWs of members in both associations we can assess any dif- 
ferences between them. It would be simplest to compare the 
EWs of likely members in the two associations using J — I 
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Table 3. Equivale nt widths an d radial velocities measured using the Nal doublet. The field number and object number refer to des- 
ignations found in lKenvon et a l. f2004). Quality flags for I and R-I are as described in Burnineham ct al. f2003). The continuum was 
measured from bands at 8151-8173 A and 8233-8255 A. The EWs were integrated over the range 8180-8202 A. Radial velocities were 
measured by cross correlating against GJ3517, and corrected to give barycentric radial velocities. 



Field 


Object 


„ / T o /"I n / 4 \ 

a(J2U0U) 


o(J2000) 


I 


o-i 


Flag 


R 


— / 




Flag 


EW 

A 


S 

A 


Vbar 

km s^"*^ 


km s^^ 


0/ 
70 


1.01 


319 


05 39 48 91 1 


-02 29 1 1 05 


15.012 


0.007 


VV 


X 


.595 


0, 


.011 


VV 


3.23 


0.11 


24.25 


0.86 


77.3 


8.02 


179 


S9 47 696 


-02 36 22 96 


15.064 


0.005 


00 


X 


.508 


0, 


.008 


00 


3.04 


0.10 


26.37 


1.13 


99.9 


4.03 


229 


05 38 23 536 


-02 41 31 66 


15.149 


0.008 


00 


I 


.686 


0, 


.013 


00 


2.31 


0.25 


27.61 


2.38 


98.8 


1.01 


253 


05 40 51 369 


-02 31 49 93 


15.232 


0.007 


00 


I 


.337 


0, 


.011 


00 


3.21 


0.30 


10.48 


4.91 


0.0 


8.02 


143 


yjo oy Ok) .'±'±0 


(19 4*^ 
-uz 00 o.to 


15.258 


0.005 


00 




.725 


0, 


.007 


00 


2.17 


0.11 


31.18 


1 00 


100 


3.01 


67 




-02 Sfl 4"^ "^S 

ou '±0.00 


15.278 


0.016 


00 


X 


.878 


0, 


.038 


00 


0.64 


0.10 


32.11 


1.32 


100.0 


2.03 


260 




09 "^^^ 9fi 8Q 
-u^ 00 ^u.oy 


15.304 


0.007 


00 


X 


.320 


0, 


.011 


00 


3.84 


0.12 


18.98 


1.12 


0.00 


4.03 


237 




09 41 Q 1 '^ 
'±j_ y.±o 


15.400 


0.006 


VV 


X 


.690 


0, 


.010 


VV 


1.32 


0.39 


33.86 


3 56 


89.5 


1.03 


60 


f)^ df\ '^f\ 1 7Q 


09 1 9 fi 1 

Oil -L U.-LO 


15.520 


0.007 


VV 


X 


.594 


0, 


.011 


VV 


4.13 


0.76 


31.87 


6 26 


84.3 


1.03 


108 




-02 10 41 98 

UiJ J-U '±-L.i'0 


15.545 


0.007 


00 


X 


.675 


0, 


.011 


00 


3.59 


0.18 


29.15 


1.93 


100.0 


8.04 


188 


AO 42 SS7 


-02 23 47 43 

^0 *^o 


15.577 


0.005 


00 


X 


4X3 


0, 


.007 


00 


3.39 


0.35 


-32.22 


2.00 


0.00 


3.01 


51 


05 38 23 070 


-02 36 49 32 

0^ OU rr^y<0^ 


15.716 


0.014 


00 


X 


.423 


0, 


.022 


00 


1.57 


0.27 


34.50 


2.51 


92.1 


2.03 


63 


do S4 389 


-02 44 9 52 


15.770 


0.007 


00 


X 


.850 


0, 


.011 


00 


3.28 


0.17 


33.15 


1.24 


100.0 


4.03 


29 


05 38 22 824 


-02 45 30 43 


15.876 


0.006 


00 


X 


3X4 


0, 


.010 


00 


3.34 


0.35 


-10.71 


7.19 


0.00 


4.03 


368 




09 "^^^ 4fi 04 


16.123 


0.006 


00 


X 


.911 


0, 


.011 


00 


2.22 


0.30 


30.38 


3.21 


99.6 


2.03 


191 


0^ 4(1 2Q 4S7 


-02 40 ^5 Q1 

rtvj 00. yj- 


16.229 


0.005 


00 


X 


.368 


0, 


.007 


00 


3.48 


0.25 


XC failed 






8.01 


333 


(T^ 4n ^^9 ^^70 


09 "^K 9'^ 4Q 


16.246 


0.005 


00 


X 


.558 


0, 


.008 


00 


2.46 


0.47 


29.98 


3 38 


99.4 


4.04 


481 


(T^ 'iS '^fi '^74 


09 47 S 99 


16.260 


0.006 


00 


X 


.527 


0, 


.011 


00 


5.35 


0.45 


11.41 


8.52 


2.0 


1.02 


237 




09 91 11 70 

Uij J- J- -1- . 1 u 


16.283 


0.007 


VV 


X 


.443 


0, 


.012 


VV 


2.32 


0.24 


-10.48 


2.10 


0.00 


4.03 


215 


uo 00 00 .oou 


09 41 Sfi 

-Uir '-tl. 00. OU 


16.472 


0.006 


00 


X 


.804 


0, 


.011 


00 


2.65 


0.39 


31.11 


4 10 


97.3 


1.01 


343 


f)^ df\ 9*^ '^SQ 


09 9^^ 97 ^^ 

Ui ^0 ^ ( . -L 


16.527 


0.005 


00 


X 


.551 


0, 


.008 


00 


5.02 


0.41 


2 29 


3 06 


0.00 


1.01 


348 




-02 28 8 1 8 


16.808 


0.008 


00 


X 


.363 


0, 


.013 


00 


3.37 


0.84 


-9.63 


17.60 


0.4 


1.02 


87 


05 39 8 076 

00 kjV O.U I U 


-02 31 32 22 


16.811 


0.008 


00 


X 


.567 


0, 


.013 


00 


3.05 


0.24 


-42.90 


5.12 


0.0 


1.02 


157 


05 39 25 246 


-02 27 48 1 5 


16.856 


0.008 


00 


X 


.550 


0, 


.013 


00 


3.64 


0.28 


28.86 


2.47 


99.8 


3.01 


480 


uo 00 00.000 


09 9^^ 1 fi*^ 

Uii .^0 -L.UO 


16.989 
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00 


X 


.777 


0, 
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00 


3.42 


0.43 


35.94 


5 78 
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77 


41 2^ fi9Q 
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L/iJ OU 00.00 


17.050 
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00 


2 


.028 
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00 
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0.85 
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110 
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Figure 8. A plot of EW(Nal) versus I-band magnitude for our 
objects with P„f,i > 80% (crosses), P^^i < 60% (open circles) and 
iMartm et al. (2004^ members of Upper Scorpius (triangles) . The 
magnitudes for the Upper Scorpius objects have been scaled for 
the distance to cr Ori. Also plotted are EW(Nal) upper limits for 
the 11 objects for which our objective measure failed to detect 
Nal. 



Figure 7. The CMD as in Figure Objects with detected 
Nal are shown as open circles. Objects with Pvel > 80% are 
filled circles, whilst circles filled with crosses are objects with 
Puel >60%. Objects for which no Nal was detected are shown 
with open triangles. The dotted line follows a NextGe n 5 Myr 
isochrone <Chabrier fc Baraffe 1 997: Ba raffe et al. 200^ . Our ex- 
pected PMS region is defined as rodward of the solid line. 



colours, since this is a fair proxy to spectral type, however 
2MASS colours are only available for the 10 brightest of our 
likely members. Instead we compare the EWs in the two 
samples using the I-band magnitudes, corrected for the dif- 
fering distances to the two associations. 

Figure|5]shows that the likel y members from our sample 
display weaker Nai EWs than the lMartm et al.' f2004) mem- 
bers of Upper Scorpius. This is not surprising in the case of 
the objects brighter than I = 17 as they will be of earlier 
spectral type than any of the' Martm et al" (20041) objects, 
and the Nal doublet gets weake r with increasing tem per- 
ature at a given value of log{g) iSchiavon et al. 19^7^ . By 
this token, we might expect the fainter objects, presumably 
with later spectral types, to display stronger Nal absorption, 
however this does not seem to be the case. Both of the likely 
members with I > 17 have EW(Nal) < 2.5 A. Inspection of 
the measured EWs for the failed detections (see Table in- 
dicates that the majority of them are consistent with values 
of EW(Nal)< 2 A. This is somewhat low er than the typi- 
cal values seen in the lMartm et aL sample. Thus we 
are presented with two options: 1) the failed detections are 
interloping giants and sub- giants with low surface gravities; 
2) many of the failed detections are member brown dwarfs 
which display weak Nal absorption. Option 1 requires many 
giant stars to be located at great distance and along our line 
of sight, yet with colours that coincide with the expected 
PMS region of C-M space. Since the line of sight towards 
(J Ori is elevated -17.3° from the galactic plane, it is hard 



to see that option 1 is feasible. Option 2, however, requires 
the low-mass objects in a Ori to display consistently weaker 
Nal absorption in the 8183, 81 95 A doublet than has been 
observed in Upper Scorpius bv lMartm et al. (20041) . 

Upper Scorpius is thought to be marginally older than 
a Ori, 5-10 Myrs versus 1-5 Myrs, so we might expect to 
see weaker Nal absorp tion in a Ori, as log(a ) will be lower 
in the younger group. iMohantv et al. (2004^ derive surface 
gravities for a number of low-mass objects in Upper Scor- 
pius, and compare them to model predictions. They find 
that log{g) changes from about 3.7 to 4.2 between 1 and 
10 Myrs for a 0.04 Mq brown dwarf (I = 18 at the dis- 
tance of a Ori), and even less of a change at lower masses. 
Given that gi ants have log (g) « 0, and display EW(Nal) 
~ 1 AlSchiav on et al. 1997^ . it is difficult to see how such a 
small change in log{g) could give rise to the factor 2 discrep- 
ancy between the EW(Nai) measured here for objects in a 
Ori and those mea sured for objects in Upper Scorpius by 
'Martm et al. (2004^. The discrepancy between the two sets 
of EWs can likely be explained without resorting to a physi - 
cal interpretation. The spectra used bv lMartfn et al. (2004^ 
are significantly lower resolution (R=900) than those used 
here, and this could lead to placing the pseudo-continuum 
over blended absorption features such as the TiO band just 
redward of the Nal doublet. Thus we suggest that option 2 
is the most likely explanation for these results. We can not 
determine if this is due to a physical difference between the 
two regions or if it has resulted from a difference in mea- 
surement techniques. 

This result implies that it will be very hard to distin- 
guish between background giants and member brown dwarfs 
based on the strength of Nal absorption. 

The low values of EW(Nal) implied for faint objects 
in our sample suggest that to avoid biasing such samples 
against finding members, observations should be made to a 
sufficient signal-to-noise to detect the Nal doublet in giant 
stars. Unfortunately, this destroys one of the main benefits 
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of using the Nal doublet for radial velocity determination in 
main sequence stars, the ability to measure it at low signal- 
to-noise due to its strength. This is offset, however, by the 
fact the gravity sensitive nature of this doublet can provide 
a second diagnostic for contaminants 

A possible weakness in use of radial velocities to de- 
fine membership of young groups is the fact that bona fide 
members can be ruled out when they are members of bi- 
nary svstems. lDolan fc Mathiexi f2001j) use the radial veloc- 
ities in concert with the presence of Lithium absorption at 
6708 A to identify 266 likely members of the A Ori young 
group. They find that 9 objects with EW(Li) consistent with 
youth are double lined binaries, while 3 appeared to be single 
line binaries. Whilst double lined binaries can be recognised 
as such, and not ruled out of membership lists, the single 
lined binaries would be ruled out in the absence of addi- 
tional evidence of their membership. If we assume a similar 
bi nary fraction and distrib ution of properties as that seen 
bv lDolan" fc Mathieu (20011) in A Ori, we can estimate that 
for our sample of 18 likely members, we would not expect 
to identify any double lined binaries, or to miss any single 
lined ones. Since we cannot use EW(Nal) to identify mem- 
bers, we cannot assess whether we have ruled out any single 
lined binaries from our membership lists. However, a signif- 
icant departure from behaviour elsewhere in the Orion OBI 
association would be required for a significant number of 
members to be ruled out due to binarity. 

We have cross-correl ated our sample against the candi- 
date membership hsts of lBeiar et al. (199911 200lh . We find 
that only one object, 3.01 480, correlates in position with 
an object, SORT 68, from their sample. However, we do not 
believe that they are the same object since there is a nearly 
6 magnitude discrepancy in their Tband brightness (1=16.99 
vs. 1=23.78). The fact that we do not share objects with the 
iBeiar et al. (1 999 2001) sample is explained by a number of 
factors. Firstly, our RI catalogue excludes the regions im- 
mediately adjacent to the a Ori multiple system, whereas 
the Bejar et al sample does not. In addition, our catalo gue 
covers a much larger Section of sky ijKenvon et al. 2004) . so 
there are many objects that Bejar et al were simply unable 
to select as candidates. Finally, we draw our sample from 
a wider region of C-M space than Bejar et al, and so we 
simply miss some of their candidates due to a lower density 
of selected objects. 



6 CONCLUSIONS 

We have carried out high resolution spectroscopy of the 
Nal doublet at 8183, 8195 A of a sample of candidate 
low-mass stellar and substellar members of the a Ori young 
group drawn from a broad region of C-M space. We have 
selected a sample of 54 objects which displayed sufficient 
signal-to-noise to detect Nal with EW(Nal) = 3 A at a 
significance of 2a, a criterion base d on observations o f 
brown dwarfs in Upper Scorpius by iMartm et al. (20041) . 
Significant (EW(Nal) > 2o-ew) Nal was detected in 38 
of the 54 sample objects, and these were cross-correlated 
against an M9V standard to obtain radial velocities, which 
were then used to calculate membership probabilities. We 
find that 13 objects are likely radial velocity members 
{P^ei > 80%) of the a Ori young group. Based on these 



probabilities, and the values measured for EW(Nal) we 
arrive at the following conclusions. 

1) Photometric selection techniques do not miss sig- 
nificant numbers of bona fide association members. 

2) At I brighter than 17 the expected PMS region of 
the CMD does not contain a significant number of contam- 
inants. 

3) Very low-mass objects in the a Ori young group appear 
to have weaker EW(Nal) than found by Martm ct al. (200|) 
for low-mass members of the Upper Scorpius OB associa- 
tion. We have no explanation for why this is so, although 
we suggest measurement effects could account for this. 

4) High resolution observations of the Nal doublet at 
8183, 8195 A offer the possibility of 2 membership diag- 
nostics for very low-mass objects from a single observation. 
Ensuring observations are made to sufficient signal-to-noise 
to detect the doublet in the spectrum of a {logg ~ 0) giant 
star will avoid biasing a sample against bona fide members 
with lower than expected EW(Nal). 
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